Budded membrane microdomains as regulators for cellular tension 
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We propose a mechanism for the control of the mechanical properties of the membrane of living 
cells that is based on the exchange of membrane area between the cell membrane and a membrane 
reservoir. The reservoir is composed of invaginated membrane microdomains which are liable to 
flatten upon increase of membrane strain, effectively controlling membrane tension. We show that 
the domain shape transition is first order, allowing for coexistence between flat and invaginated 
domains. During coexistence, the membrane tension is controlled by the domains elasticity and by 
the kinetics of the shape transition. We show that the tension of the plasma membrane of living cells 
is inherently transient and dynamical, and that valuable insights into the organization of the cell 
membrane can be obtained by studying the variation of the cell membrane tension upon mechanical 
perturbation. 
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Endocytosis, exocytosis, cell motility and many other 
crucial cellular processes are known to be influenced by 
the tension of the cell membrane 0. The force needed 
to pull a tubular membrane tether from the cell with an 
optical trap is often used as a probe of the mechanical 
state of the membrane^ . The mechanical response of liv- 
ing cells to such perturbation involves cytoskeleton defor- 
mation, the breaking of membrane-cytoskeleton bondsQ 
and changes in membrane morphology. The level of mem- 
brane tension in cells is thought to primarily reflects cy- 
toskeleton anchoring to the membrane J2jj. Nevertheless, 
there is experimental evidence that Surface Area Regu- 
lation (SAR) occurs Q| and is able to buffer the increase 
of tension upon mechanical perturbation. SAR could be 
achieved by the transfer of membrane area from a reser- 
voir to the plasma membrane, during which the level of 
membrane tension is at least partly controlled by the 
dynamical response of the reservoir. One possible man- 
ifestation of the exchange of membrane with a reservoir 
is the presence of a plateau in the force-extension curve 
of tether extraction experiments, as observed in The 
plateau ends when the putative membrane reservoir is 
emptied. In this particular experiment the reservoir is 
able to react quite fast to the perturbation (< 0.1s), 
which could indicate i) that the reservoir is permanently 
connected to the plasma membrane and ii) that the reg- 
ulation is purely physical and does not involve biological 
signaling. 

Our goal in this paper is to study the influence of 
membrane morphology, and in particular the flattening 
of invaginated "bud- like" domains, on the mechanical re- 
sponse of the membrane. Such domains include caveolae 
that form fl shaped invaginations reactive to mem- 
brane stress0, Striking experiments, showing Cave- 
olae flattening under tension 9] , support the idea that 
the delivery of invaginated membrane area to the plasma 



membrane is controlled by tension [lfl|. The shape of 
membrane domains in artificial membrane systems, such 
as giant vesicles, is also known to be dependent on the 
membrane tension [lj. Indeed, domain budding has been 
observed upon the decrease of membrane stress [l^. 
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FIG. 1: Sketch depicting the mechanical response of a mem- 
brane to the pulling of a membrane tether. Pulling an area At 
out of the cell either increases the membrane tension 7 (the 
length of the spring), or triggers membrane exchange with a 
reservoir (flattening of invaginated membrane domains) 

Our model is summarised in Fig^ Here the cytoskele- 
ton is included as linear elastic element (spring) mechan- 
ically coupled to the cell membrane. Although crude, 
this approximation is consistent with the experimental 
results of (see 0). Furthermore, it allows us iso- 
late and study the response of the reservoir, which can 
be thought of as one element in a chain of response 
that transmits stress to the cytoskeleton. Cytoskeleton 
elasticity and membrane cytoskeleton anchoring are re- 
flected by an effective stretching coefficient K s of order 
Within this framework, the variation 
of the cell tension with the tether area At reads: 



7 = 70 + K 8 (At + Ares ~ A r °J s ^ 



(1) 



where 70 and A r °e\ are the cell tension and reservoir area 
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at rest (without tether, At = 0). If membrane area is 
delivered to the plasma membrane, as during exocytosis, 
the "tether" area is negative. One can see that the mem- 
brane tension 7 can be maintained constant upon tether 
pulling only if the decrease of reservoir area matches the 
increase of tether area. 

The membrane domains constituting the reservoir are 
described within a very general theoretical framework. It 
is based on the domain bending rigidity, that disfavors 
the budded state, and its composition difference with the 
rest of the membrane, that promotes domain invagina- 
tion in order to reduce the length of the domain periph- 
ery (FigGJl. Domains are treated as spherical caps of 
fixed area S and adjustable curvature C . Their shape is 
uniquely characterized by the parameter f3 (FigEl an< i 
Eq.(0l), equal to unity for a fully budded domain (a 
sphere), and that vanishes for a flat domain. The energy 
/ of a domain contains a surface tension-independent 
part This itself involves a term arising from the 

line tension a , proportional to the length of the cap 
edge (neck), and a term giving the bending energy of 
the cap, proportional to the bending rigidity K and the 
squared curvature [L4|. Simple geometry then gives 

U tK \0\ = V^SayJl - [3 + 8ttk/3 [3 = SC 2 /16ir (2) 

If the domain size exceeds a critical value S c = 7t(4k/(t) 2 , 
an invaginated sphere (/3 = 1) has a lower e nerg y than a 
flat domain (/3 = 0) and budding is expected |l5|. We as- 
sume that S > S c in what follows. Typical values of the 
parameters a ~ ksT/nm and k = 20/csT correspond to 
a critical size S c = (lOOnm) 2 , similar to the size of caveo- 
lae. In practice invaginated domains remain attached to 
the mother membrane by a small neck. Here, this is con- 
trolled phenomenologically by assigning the value otbud 
(3> 1) to the ratio of invagination radius R to neck size, 
giving an upper bound fibud = 1 — (c*6«d/2) 2 to the shape 
parameter. 




FIG. 2: Left: A sketch of a flat and an invaginated domain, 
above our idealization for the domain shape; a spherical cap of 
area S, radius R (curvature C = 2/R). The shape is defined 
by a single parameter < f3 — SC 2 /16tt < 1. Right: energy 
of a domain for increasing value of the membrane tension. 
Low tensions (red) favor budded shape (/3 = 1), and high 
tensions (blue) favor flat shapes (/3 = 0). At coexistence 
(green), flat and budded domains have the same energy, and 
the tension is given by Eq.l|3J The budded state remains a 
local minimum for high membrane tension and there exists 
an energy barrier A/ to domain flattening. 



Including the membrane tension 7, the domain energy 
reads: f[f3] — f a , K [0\ +jSf3. Increasing membrane ten- 
sion increases the energy of curved states {3 > and pro- 
motes the flat state (see FigEJ. The flat states eventually 
becomes stable for a critical tension 7*: 

S7* = -(/<.,«[!] - = 2v^5a - Snn (3) 

As can be seen in FigEl the budded and flat domain 
shapes are separated by an energy barrier for interme- 
diate tension. The existence of this barrier is crucial to 
their function as tension regulators, as it allows the co- 
existence of flat and invaginated domains. Variation of 
the membrane strain (the tether length) may occur with- 
out change of tension by adjusting the fraction of budded 
domains. 

The energy scales in this system (with S = O.l^tm 2 ) 
are 

a = VttSct ~ R = 87tk ~ 500fcsT ; 7 = (4) 

The energy of surface tension competes with the line and 
bending energies for 7 ~ 500fcsT, or 7 = 2 10~ 5 J/ra 2 
(corresponding to tether forces of order WpN). This 
is precisely in the range of mechanical tension recorded 
for cellular membranes |2(, which is very encouraging for 
the biological relevance of our model. One notes that 
the energy scale is very large compared to the thermal 
energy ksT, or to the energy of any "active tempera- 
ture" present in biological sy stems [HI . This has two 
important physical consequences: (i) the shape transi- 
tion of a domain is very discontinuous, a domain snaps 
open rather than continuously flattening upon tension 
increase, and (ii) the budding and flattening transitions 
should actually occur at different tensions, for which the 
respective energy barriers are of order ksT. In biolog- 
ical systems, the "temperature" T might be seen as a 
parameter reflecting cellular activity, such as the poly- 
merization of the actin cortex near the membrane, and 
the activity of membrane pumps. For simple cells such 
as Red Blood Cells, it is typically a few times the ther- 
modynamic temperature 

The bottleneck for the shape transition is the maxi- 
mum of energy, which corresponds to a shape parameter 
Pmax — 1 — (^/(ft + l)) 2 ■ The budding and flattening 
tensions (7 1 - 1 -' and respectively) at which the corre- 
sponding energy barrier vanishes are: 

7« = ct-k<7*=2ct-k< 7 (0) = 2a bud a - R (5) 

where oibud characterize the size of the invagination neck 
(see above). 

We investigate the tension regulation performed by a 
collection of M domains, of total area NSflbud (where 
Pbud = 1 — (afcud/2) 2 , is the largest value of shape param- 
eter consistent with the existence of a finite size neck). 

When flat and budded domains coexist, a fraction e 
of domains are invaginated, and the reservoir area is 
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Ares = SNeflbud- The total membrane energy, includ- 
ing the contribution f aiK of each of the J\f membrane 
domains (Eq.(J2J) and the total work done against mem- 
brane tension can be written: 

T = M (efa, K \j3bud\ + (1 - e)f*,*[°]) + J dA^A] (6) 

Optimizing the energy for the fraction of invaginated 
domains dT '/de = leads directly to the regulation of 
membrane tension, when flat and budded domains coex- 
ist (0 < e < 1). Substituting in Eq.© the expression 
for the surface tension Eq.QJ (with ftbud ^ 1), we find 
that the tension is set to the value 7* of Eq.©, which 
depends on the characteristics of the membrane reservoir 
(cr, k, and S), but not on the tether area. Regulation is 
achieved by adjusting the fraction of budded domains to: 



(70 + K S A T /S) - r 
K S M 



K s = K S S 2 



(7) 



where eo is the fraction of budded domains corresponding 
to the tension at rest 70 , and where a normalized stretch- 
ing coefficient K s (^ O.Ict), with dimension of energy, is 
introduced for convenience. 

If the reservoir is given time to equilibrate, regulation 
starts for a level of perturbation corresponding to a tether 
area Aj^ (at which all domains are budded, e* = 1), and 
ends at a tether area A^ (at which all domains are flat, 
e* = 0). The tension of the cell membrane is then set 
to the values 7*, for any perturbation within the range 
A£ < At < A^! . If the perturbation is very fast, one 
expect a large difference between the regulated tension 
upon tether pulling and tether retraction, in agreement 
with Eq.JSJ). 

To obtain the full kinetic response of the membrane to 
strain, we describe the transition as a classical Kramers' 
process [l^. were the transition time between two states 
is exponential with the energy barrier A/ that has to be 
overcome in the process: r = tq exp[A//fcsT]. Here, To 
is the characteristic fluctuation time of the domain shape, 
and T may be an effective temperature resulting from the 
large biochemical activity near the cell membrane[lfij|. 
The transition time is very much dependent upon the 
membrane tension. Assuming that the transition of a 
single domain occurs with negligible change of tension 
(this implies M ^S> 1), the transition is fully described by 
the energy /[/?] = / CT>K [/?] + 7/3, with 7 given by Eq.(HJ. 
Assuming for simplicity that both domain flattening and 
budding involve the same fluctuation time To, the kinetic 
evolution of the fraction e is given by 



de 
T °di 



f m ax-fW bud ] 



+ (l-e)e- 



-/[ii] 



(8) 



where the maximum of energy f max corresponds to the 
least favorable domain shape (3 m ax = 1 — a 2 /(k + 7) 2 



In order to mimic the tether pulling experiment, 
where the tether is typically extracted at constant speed 
(4/Ltm/s) in |3| , we consider the reservoir response to 
a perturbation applied with a given rate At- At = 
A^+ATt. If the perturbation is applied slowly (Attq <C 
Ares), the reservoir has time to equilibrate (de/dt = 0), 
and the fraction e T of budded domain is found from 



Eq.© to be given by e T /(l - e\ 



M[0]-f^ d )/(k B T) 

with/[0] — fb u d — 2a — (7(e) + R). The fraction eT is 
the equivalent of the equilibrium fraction e* (Eq.Q)), 
that takes thermal fluctuations into account (e^ ~ e* 
if d 3> ksT). Thermal fluctuations smoothen the transi- 
tion between budded and flat domains by allowing states 
of non-minimal energy to be populated. As a conse- 
quence, the tension is not perfectly constant during the 
transition, and the slope at mid plateau is of order 
dj*/dA T \ plat ^ Ak B T/(SAres)- If, on the other hand, 
the perturbation is applied very fast, the shape transition 
requires small energy barriers, which means high tension 
for bud flattening, and low tension for domain budding, 
respectively 7^ and 7W give by Eq.JSJ). 
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FIG. 3: Sketch of the change in membrane morphology upon 
tether extraction. At low strain (1) all membrane domains 
are budded, and the membrane tension 7 increases linearly 
with the membrane area At- For large perturbation (3) all 
domains are flat and a similar linear increase is observed. In 
between (2) the membrane tension is maintained at a plateau 
value while flat and budded domains coexist. The equilib- 
rium reservoir response (7*, dashed gray line) corresponds 
to a quasi-static perturbation. The dynamical responses at 
constant rate At is shown upon extraction (red) and retrac- 
tion (blue) for various rates. The hysteresis, which increases 
with the perturbation rate, illustrates the kinetic nature of 
the domain shape transformation 

The physical mechanism at the origin of tension reg- 
ulation and the membrane hysteretic response to tether 
extraction and retraction are shown in Fig|3] To obtain 
an analytical expression of the plateau height with the 
perturbation rate, we approximate that the tension is al- 
most constant during regulation (d'y/dt ~ 0) so that the 
energy barrier is of order A//(fcsT) ~ log A res / {Attq) . 
The plateau tension upon increase and decrease of the 
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perturbation are then respectively given by 71^0 



v(0). 



3/2 



\J crkgT log and 70^1 = T ' + \fak B T log, 



7^-2a b „ d 

with log = log[yl res /(2^4TTo)]- As expected for acti- 
vated processes, the dependence of the tension at tran- 
sition with the rate of perturbation At is logarithmic. 
The same is true for the slope of the plateau, which 
can be estimated by identifying the plateau inflexion 
point. The condition d 2 ^/dt 2 = imposes A//(&bT) ~ 
e/e, corresponding to a plateau slope d^ji^o/ dAr\ plat — 

y / '2kBT(jabud 3 1 log/*4res- As the effective membrane 
temperature may be inferred independently from mea- 
surement of membrane fluctuations [l6|, the study of both 
the height and slope of the tension plateau gives valuable 
information of the kinetics of area transfer between the 
reservoir and the plasma membrane. 

The difference between the quasi-static and dynamic 
plateau can be of order 1CP 4 J/m 2 and correspond to a 
difference in force of order 10p./V. This is precisely the 
scale of the forces measured upon tether extraction in 
This possibility thus exists that the initial increase 
of tension observed prior to the plateau in || is of purely 
kinetic origin and originate from the slow response of an 
already partially unfolded reservoir at rest. In this case, 
the membrane tension of the cell at rest would be 7* of 
Eq.©, fully controlled by the mechanical properties of 
the membrane domains forming the reservoir 

In summary, we have derived the mechanical reactiv- 
ity of a cell membrane, in contact with a reservoir com- 
posed of invaginated membrane domains liable to flat- 
ten under strain. The flattening transition is first or- 
der, which means that the invaginations snap open above 
a critical strain rather than continuously flattening, an 
observation consistent with experimental evidences from 
the membrane invagination caveolae|lfjj. As a conse- 
quence, the cell mechanical response shows a plateau 
during the transition, corresponding to the coexistence 
of flat and invaginated domains. This study provides the 
basis for a mechanical regulation of the tension of the cell 
membrane. The mechanism at the origin of this regula- 
tion is of purely physical origin, consistent with the fast 
timescale (< O.lsec) of the observed cellular response. In 
practice, a ring of specialized membrane proteins such as 
dynamin^^ is often present at the neck of membrane in- 
vagination. These proteins most probably influence the 
domain line energy, and might even dominate the energy 
required to flatten the domain. The regulation of mem- 
brane tension by membrane invaginations rely on the ex- 
istence of two well defined domain shapes, separated by 
an energy barrier. If anything, neck proteins can only 
increase the energy barrier to flattening, thereby rein- 
forcing tension regulation. 

This work also opens the possibility of a new quanti- 
tative "force spectroscopy" of the cell membrane. One 



could thereby obtain structural information on the mem- 
brane organization, in much the same way information on 
a protein structure can be gathered from force measure- 
ment upon protein unfolding As a first step, one 
may identify the fairly regular oscillation of the force 
during regulation in ; 5j, to the flattening of single do- 
mains. Preliminary analvsis|l3j hints at domains of area 
S ~ (400nm) 2 . 
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